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By combining our theoretical calculation and recently measured electric quadrupole hyperfine struc- 
ture constant of the 3d ^1)5/2 state in the singly ionized ^^Ca, we determine its nuclear quadrupole 
moment to one percent accuracy. The obtained result, —0.0444(6)6, is about ten percent improve- 
ment over the considered standard value. We have employed the relativistic coupled-cluster theory 
at single and double excitations level to calculate the electronic wave functions. The accuracy of 
these wave functions are estimated by comparing our calculated magnetic dipole hyperfine constants 
with their corresponding available experimental results of many low-lying states. We also present 
hyperfine structure constants for other higher excited states where experimental results are not 
reported. Role of the Breit interaction has been investigated in these properties. 

PACS numbers: 21.10.Ky,31.10.+z,31.30.Gs,32.10.Fn 
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I. INTRODUCTION 

Advanced modern techniques of laser cooling and trap- 
ping have enabled us to carry out precision measurements 
of hyperfine structure constants in atomic systems P, ■ 
Theoretical studies of these quantities require accurate 
many-body methods, inclusion of relativistic effects and 
knowledge of nuclear moments [1, 0, |^, @ . Precise mea- 
surements of nuclear moments are difficult, especially the 
quadrupole and octupole moments. The prominent ex- 
amples of techniques to measure them are NMR, atomic 
beams, optical pumping, recoil methods etc. How- 
ever, their absolute results are of great interest for the 
nuclear ph ysic ists to be able to test different nuclear 
models [3, [l3|. The quadrupole moment of the stable 
isotope ''^Ca is also of particular interest in the evalua- 
tion of the nuclear ma gnet ic resonance measurements in 
the biological systems [111 [l^ . Investigating properties 
of nuclei in the region of magic numbers are challeng- 
ing because the valence nucleons can be strongly affected 
by the close shell configuration. ''^Ca has mass number 
between the double magic numbers ''°Ca and ''^Ca sug- 
gesting peculiar nuclear charge distribution. It is plau- 
sible to obtain correct electronic wave functions, hence 
their properties, in the single valence atomic systems 
using the recently developed highly potential methods 
like all order relativistic coupled-cluster (RCC) theory 
SlilQIil- Singly ionized ^^Ga. {^^Ga + ) is of medium 
size system and single and double excitations approxima- 
tion RCC approach (CCSD method) is capable enough 
in such cases 0, [3, to account electron correlation 
effects accurately. It has been observed that electron cor- 
relation effects exhibit spectacular behavior in the stud- 
ies of the magnetic dipole hyperfine structure constant in 
the ^Dqi2 states of the singly ionized alkaline-earth metal 
atoms. The CCSD method with contributions from lead- 



ing order triple excitations (CCSD(T) method) is able to 
consider them sufficiently for which it can produce results 
matching with the precisely measured values [l^, [l^, [13] . 
In this work, we have used this method to evaluate the 
electronic matrix elements due to the hyperfine interac- 
tion operators. 

It would also be imperative to mention here that '*'^Ca"'" 
is an interesting candidate for quantum computation [l^ 
and optical frequency standard [l^ . Accurate values 
of hyperfine structure constants are useful in estimating 
shifts in the energy levels due to the stray electromagnetic 
fields during the experimental set-up [2l|,[22]. Again, the- 
oretical estimations of these quantities are used to test 
the correct behavior of wave functions in the nuclear re- 
gion [H,[23|. Nuclear magnetic moment of ''^Ca is known 
within sub-one percent accuracy [5^ , but its reported nu- 
clear quadrupole moment values vary from —0.0408(8)6 
to -0.065(20)6 [H [13, [H, [H, [33. Using the accurate 
matrix element of the electric quadrupole hyperfine in- 
teraction operator and precisely measured [31] electric 
quadrupole hyperfine structure constant of the id 
state in ^^Ca+, we determine its nuclear quadrupole mo- 
ment. Accuracy of the theoretical calculations are esti- 
mated from the analysis of the correct behavior of the 
wave functions which are able to reproduce the magnetic 
dipole hyperfine structure constants in a few low-lying 
states where experimental results are available. We also 
present hyperfine structure constants of other higher ex- 
cited states for further testimony of our results using fu- 
ture experiments. 



II. THEORY AND METHOD OF 
CALCULATIONS 
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The detail theory about the hyperfine structures are 
given earlier in a classic paper by C. Swartz [^. Here, we 
have mentioned only these formulae in explicit form. The 
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relativistic hyperfine interaction Hamiltonian is given by 
if,,/, = Vm«.T« (2.1) 



where M*^') and T^'^ are spherical tensor operators of 
rank /. In the first-order perturbation theory, hyperfine 
interaction energy Wf of the hyperfine state \F; I, J) of 
angular momentum F = I -\- J with / and J being the 
nuclear spin and electronic angular momentum of the 
associated fine structure state \ J,Mj), respectively, after 
neglecting terms beyond I = 2 is given by 



B 



hfs- 



1) - 2/(7+ 1)J(J+ 1) 



2/(2/- 1)2J(2J- 1) 



(2.2) 

where K = F{F + 1) - I{I + 1) - J(J + 1), Ahjs is the 
magnetic dipole structure constant for I = 1 and /J^/s 
is the electric quadrupole structure constant for I ~ 2. 
These constants are defined as 



with the angular momentum selection rule 7^(1 f, I, h) = 1 
when If + I + li =evcn for the orbital angular momenta 
If and li, otherwise zero. 

The primary objective of this work is to calculate the 
above electronic matrix elements of the hyperfine inter- 
action operators. It is obvious from the single particle 
expressions that these matrix elements have strong over- 
lap with the nucleus. As a starting point, we consider 
kinetically balanced Gaussian type of orbitals (GTOs) 
which seem to be an ideal choice for obtaining correct 
behavior of wave functions in the nuclear region [33 , [ssj 
to calculate mean- field wave functions j^o) of the closed- 
shell configuration in the Dirac(Hartree)-Fock (DF) ap- 
proach. To calculate the atomic state function (ASF) of 
single valence with closed-shell configurations, we express 
it in the RCC ansatz as 



(2.9) 



A 



hfs — ^J'NgI 



(J||T«||J) 



^J(J+1)(2J+1)' 



and 



B 



hfs 



8J(2J- 1) 



(2J+ l)(2J + 2)(2J + 3) 



1/2 



where is the new reference state which is defined as 

$„) = aJ|$o) and will give the DF result for the above 

(2.3) open-shell configuration. In the above expression, we call 

T and Sv as the closed and open shell core and core 

with valence electron excitation operators, respectively, 

which in the second quantization notation in the CCSD 

/ON approximation are given by 

(J||t(2)||J},(2.4J'^ ^ 



where we have used atomic unit (au). In the above ex- 
pressions, fiN is the Bohr magneton and we use gj = 
[^] with /i/ and I{— 7/2) are the nuclear magnetic 
dipole moment and spin, respectively, as —0.37646943 
from the measurement [25| to evaluate A^fs- Since nu- 
clear quadrupole moment, Qnuc, is not known accurately, 
hence we calculate 



T^T,+T2^Y1 + \Y1 «a6aaiS?,(2.10) 

a,p ab,pq 

Sv = Sly + S2v = a^avsl + ^ X! '^p^q'^bavsl%^ .11) 



b,pq 



B 



hfs 



Q7} 



8J(2J- 1) 



(2J+ l)(2J-|-2)(2J + 3) 



1/2 



(J||t(2)||J).(2.5) 



The reduced matrix elements of the electronic spherical 
operators, T^'-* = X^^''^ ii^ terms of single orbitals are 
given by 



(^/||iW||«.) = -(^/ + ^.)(-«/l|CW||«.) 

' ^^{PfQ^±Qf31 



(2.6) 



and 

{nf\\t('^\\n,) = -{nf\\C^'^\\^,) 



dr 



(PfPr+QfQr) 



, (2.7) 



where the (a, 6, c), {p,q,r..) and [v) subscripts of the 
second quantized operators represent core, particle (vir- 
tual) and valence orbitals, respectively. The t and Sy 
coefficients are the corresponding excitation amplitudes 
which are determined using the following equations 

($^|{i^^}|$o) =0 (2.12) 

with the superscript L(~ 1, 2) representing the single and 
double excited states from the corresponding reference 
states and the wide-hat symbol denotes the linked terms. 
A/Jy is the corresponding valence electron affinity energy 
which is evaluated by 



where i and / represent initial and final orbitals, respec- 
tively. The reduced matrix elements of the spherical ten- 
sors (C«) are given by 



AEy = (^y\{HNeT}{l + S4\^y) 



(2.14) 



{^f\\C'^'^\\'^^) = (-l)^/+i/y(2j/ + l)(2j, + l) 

Trilf,l,k) (2.8) 



In our CCSD(T) approach, we consider effects of the 
leading order triple excitations through AEy by con- 
structing triple excitation operator 



Jf 
1/2 



-1/2 



qpqr _ 



HnT2 + //at 5*1,2 
Cfc + — £(7 — Ct 



(2.15) 
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and contracting S^l^ with the Hamiltonian to get con- 
tributions to the corresponding AE^, where is the DF 
energy of the electron in the i*^ orbital. 

We consider the Dirac-Coulomb-Breit Hamiltonian in 
the above equations that is given by 



TABLE I: Results of Ahfs, Bhfs/Qnuc and Bhfs of many 
states in ''^Ca^. 



1 



F = ca ■ p + (/? - ly + Vr,uc{r) + 

1 \ ai-Q.2 _ {3.1 ■ ri2)(a2 • ^2) ~ 
2 



+ 





Ahfs 


-^hfs / Qnuc 


Bhfs 


Calc 


Expt 


Calc 


Expt 


(MHz) 


(MHz) 


(MHz6-i) 


(MHz) 



^12 



r 



12 



%16) 



where c is the velocity of light, a and P are the Dirac 
matrices and Vnuc{r) is the nuclear potential. We solve 
the wave functions due to the above Hamiltonian in the 
DF and RCC methods self-consistently with the tolerance 
size below 10^^ to obtain precise results. 

We evaluate expectation values due to the hyperfine 
interaction operators using our RCC method by 



<0> = 



<a>,|{i + st}0{i 



4s "51/2 
5s 

4p 'Pl/2 
5p 'Pl/2 

4p 'P3/2 



<$„|{1 + 5'I, + 4JO{1 + 5i„ 



{1 + si + SlJNoil + Si, + S2v} 



17) 



where O is any of the operator, O = [e^ Oe^) and 
iVo = e^^e^. Generally, both O and in the RCC 
approach are non-terminative series. However, we use a 
special trick to obtain their almost all the leading order 
contributions using the Wick's generalized theorem 34 1. 
In this procedure, we evaluate first effective one-body, 
two body terms etc. step by step and at the end sand- 
wich them between the Sv and its conjugate operators. 
This procedure has already been demonstrated in our 
earlier works [l^, IH, |3l, [Sg . We also explicitly present 
contributions from the normalization factors evaluating 
them in the following way 



Norm = (^'„|0|*„){- 



1 



-1}, 



where iV„ = {1 + S^ + 4j^o{l + ^1. + ^2.}- 



5p 'P3/2 



■^3/2 



4d 2_D. 



3/2 



3d ^Dr^i2 



-806.387(1.0)" 
-805.348*' 
-8ir 
-794.7'* 

-234.0(2.0)" 

-145.422(340)" 
-143.068'' 
-148"= 
-144.8'' 

-49.651(521)" 

-30.370(415)" 
-30.498'' 
-30.9'= 
-29.3'* 

-10.335(112)" 

-47.282(289)" 

-47.824'' 

-52"= 

-49.4'* 

-47.27^' 

-9.474(042)" 

-3.622(327)" 

-3.552'' 

-5.2'= 

-4.2'* 

-4.84^ 



(2.18) 4d -2.971(035)° 



-797.5(2.4)'= 

-805(2)-'' 

-817(15)" 



-158.0(3.3)'= 
-145.5(1.0/ 
-142(8)'' 
145.4(0.1)' 



-29.7(1.6)"= 

-31.9(2)^ 

-31.0(2)' 



-48.3(1.6)'' 
-47.3(2)' 



-3.8(6)' 
-3.8931(2)*= 



-6.7(1.4)^ 
-6.9(1.7)' 



-3.7(1.9)' 



151.246(752)" 
151.798'' 

155"= 



50.904(571)" 

67.280(688)" 
68.067'' 

68'= 

72.06^ 

17.442(134)" 



95.472(1.276)" -3.9(6.0)' 

100.208'' -4.241(4)*= 
97'= 

102.45^ 
24.769(213)" 



"This work; ''^8 



I; "HI; "HI; 'gj; 



III. RESULTS AND DISCUSSIONS 

Earlier, we have studied behavior of the electron corre- 
lation effects in the magnetic dipole hyperfine structure 
constants in the considered system using the CCSD(T) 
method and GTOs for few low-lying states [s^ and later 
only in the 3d state fisll. Due to the limitations 

with the available computational resources at that time, 
we had restricted our calculations by considering only up 
to f-symmetry (orbital quantum number I = 3) orbitals. 
Again, we had used basis along with g-symmetry orbitals 
in the latter to find out the peculiar behavior of the core- 
polarization effects in the ^£>5/2 states of the alkaline 



earth-metal ions, however detail investigations in the ac- 
curacy of the wave functions for other low-lying states 
were not carried out in that work. In the present work, 
we use larger basis functions up to g-symmetry that pro- 
duces the other properties [13] including electron affinity 
energies of many low-lying states matching with the ex- 
perimental results. 

In Table [H we present our Ahfs and Bhfs/Qnuc (or 
Bhfs) results of many states along with the available 
theoretical values and experimental measurements. Es- 
timated errors in our calculations are given inside the 
parenthesis. Considered sources of these errors are ba- 
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TABLE II: Breit contributions (Abt) to the 
Bhfs/Qnuc in MHz and MIIz6~^, respectively. 



A 



hfs 



and 



State 










Ahfs 




-Bfifs/ Qnuc 



4s 'Sy2 

5s 

5p'Pi 



4p 
5p 
3d 
4d 
3d 
4d 



/2 
P3/2 
P3/2 



3/2 
'3/2 
'5/2 
'5/2 



-0.697 
-0.214 
-0.115 
-0.063 
0.048 
0.014 
-0.102 
0.135 
-0.102 
-0.035 



-0.242 
-0.054 

0.150 
-0.231 

0.385 
-0.298 



sically two folds: (a) numerical calculations with limita- 
tions over finite size basis functions within g-symmetry 
and (b) approximations at the level of excitations in the 
RCC approach. In the first case, we have tested our DF 
results for a set of basis functions to achieve consistent 
results and possible discrepancies in these results are as- 
sumed as one of the sources of errors. Second, we approx- 
imate our level of excitations at the singles and doubles as 
Ca"*" is assumed as not too heavy system. However, equa- 
tions to determine these amplitudes, in principle, should 
couple with the higher excitations for accurate calcula- 
tions. Although our leading order triple excitations take 
care of most of these contributions, wc estimate the con- 
tributions from the higher excitations by studying differ- 
ences of results between the CCSD and CCSD(T) meth- 
ods and evaluating lower order diagrams that may arise 
through the neglected triple excitations in the RCC ap- 
proach. This procedure may not be sufficient enough, 
however it explains their importance qualitatively. We 
have scaled these contributions to account as the upper 
values of the second source of errors. 

It is obvious that our A^fs results match quite well 
within the uncertainties of available experimental values 
giving an indication that our calculated wave functions 
are accurate enough in the nuclear region. There are also 
experimental results available for Bhfs in many cases, but 
none of them are accurate enough except the recent mea- 
sured value in the 3d '^£'5/2 state (31j . It is obvious from 
Table U that there are quite well agreement between dif- 
ferent calculations of Bhfs/Qnuc results at least in the 4p 
^^3/2 and 3d ^0^/2 states, but there comes large discrep- 
ancies between the results in the 3d state. Due to 
consistency among the theoretical calculations in the 4p 
and 3d ^l?3/2 states, it would indeed be appropriate 
to combine these results with the measured Bhfs values 
of their corresponding states to determine less accurately 
known Qnuc in ^^Ca. In contrast, Bhfs of the 3d ^1)5/2 
state is measured quite precisely and hence it is necessary 
to use its calculated Bhfs/Qnuc result to combine with 



its measured Bhfs value to determine Qnuc in the con- 
sidered system. Therefore, we would like to investigate 
possible reasons of the discrepancies among the theoret- 
ical methods which are employed in these calculations. 
First, we investigate the role of the Breit interaction from 
which we can realize the effect of higher relativistic ef- 
fects in the above properties which was not considered 
in the previous works, then we proceed with describing 
differences in the inclusion of various electron correlation 
effects through the employed theoretical methods. 

For high accuracy calculations, it may be necessary to 
find out contributions from higher order relativistic cor- 
rections than the Coulomb interaction; occurs through 
the exchange of longitudinal photons. The next impor- 
tant contribution from frequency independent Breit in- 
teraction due to the transverse photon [i^ can be as- 
sumed as a bench mark test to estimate how big the 
neglected relativistic effects would be. In Table [TTl we 
present contributions from the Breit interaction in the 
Ahfs and Bhfs/Qnuc calculatious for different states. As 
seen, these contributions in the considered system are 
not large and it is larger in the ground state than in the 
excited states. It seems from this study that higher or- 
der relativistic corrections like bound state QED effects 
are not important for high precision calculations in the 
present system and the results mostly depend on the elec- 
tron correlation effects. 

We now discuss the differences between various works 
that account correlation effects at various level of orders. 
As mentioned earlier, our previous and present works are 
carried out with all order RCC method, but the main dif- 
ferences in the results are due to inclusion of orbitals from 
g-symmetr y a nd Breit interaction in this work. Both 
Yu et al. [381 and Martensson et al. [s^ have carried 
out their calculations using semi-empirical feature many- 
body methods. They contain all order core-polarization 
effects, however other correlation effects like Bruckner 
pair correlation effects [13] are taken up to certain orders. 
In the work of Yu et al., they have restricted orbitals in 
the evaluation of the hyperfine structure constants for 
individual state by selecting maximum contributing an- 
gular momentum configurations. However, contributions 
from all these orbitals are intrinsically accounted through 
the coupled equations in the RCC method. In Tables Hill 
and IIVI we present contributions from individual RCC 
terms to our Ahfs and Bhfs/Qnuc results, respectively. 
As seen our DF results match well with the Martensson 
et al., but there are differences between all order core- 
polarization contributions between their work and ours. 
Our all order core-polarization effects are associated with 
the S2V RCC operator [3^ [36| . The reason of discrepan- 
cies could be due to the fact that all correlation effects 
are coupled in the RCC method in contrast to the above 
semi-empirical features. There is also one non-relativistic 
theory with relativistic corrections under random phase 
approximation (RPA) approach has been employed (40| 
for calculating these hyperfine structure constants. Re- 
cently, another calculation has been carried out using 
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TABLE III: RCC contributions to the A^fs calculations. 
RCC terms 4s "gi/a 5s ^gya 4p 5p ^Pya 4p ^Pg/a 5p "P3/2 3d 4d ''0^2 3d ^1)5/2 4d "£)5~ 



(DF) 


-587.902 


-181.120 


-101.559 


-36.396 


-19.669 


-7.056 


-33.409 


-8.104 


-14.235 


-3.455 


0-0 


-1.626 


0.111 


1.076 


0.370 


0.215 


0.080 


-0.604 


-0.079 


-0.283 


-0.041 


OSi^ + cc 


-103.321 


-21.064 


-21.032 


-5.887 


-4.089 


-1.159 


-8.373 


-0.244 


-3.554 


-0.100 


052^ + cc 


-102.224 


-28.843 


-19.997 


-6.345 


-5.386 


-1.812 


-3.120 


-0.926 


15.803 


1.138 




-4.527 


-0.613 


-1.113 


-0.245 


-0.217 


-0.049 


-0.571 


-0.031 


-0.241 


-0.013 


SlOS2v+CC 


-7.338 


-0.873 


-1.689 


-0.313 


-0.374 


-0.048 


-0.260 


0.115 


0.932 


-0.112 


Sl^OS2v + cc 


-9.615 


-3.117 


-1.108 


-0.355 


-0.837 


-0.348 


-2.017 


-0.384 


-2.125 


-0.443 


Norm 


8.466 


1.519 


1.107 


0.273 


0.232 


0.057 


1.072 


0.179 


0.081 


0.055 



TABLE IV; RCC contributions to the B^fs/Qnuc calculations. 
RCC terms 4p "P3/2 5p ^P3/2 3d ''D3/2 4d "P>3/2 3d 4d 



(DF) 


96.976 


34.789 


55.354 


13.354 


78.466 


18.939 


0-0 


-0.869 


-0.307 


3.208 


0.493 


4.456 


0.680 


OSiv + cc 


20.220 


5.736 


15.275 


0.306 


21.509 


0.451 


OS2V + cc 


32.885 


10.217 


-5.066 


2.836 


-6.843 


4.067 




1.076 


0.243 


1.155 


0.060 


1.619 


0.084 


SlOSiv + cc 


1.770 


0.286 


-1.586 


0.456 


-2.225 


0.632 


SL0S2v + cc 


0.341 


0.220 


0.621 


0.308 


0.844 


0.432 


Norm 


-1.153 


-0.280 


-1.681 


-0.371 


-2.354 


-0.516 



TABLE V: Comparison of Q„uc 


values (in b) from various 


works. 




Qnuc 


Reference 


-0.0444(6) 


This work 


-0.040(8) 


[10] 


-0.0408(8) 


[26] 


-0.043(9) 


[30] 


-0.065(20) 


[28] 


-0.049(5) 


[29, 47] 


-0.044 


[38] 


-0.062(12) 


[48] 



TABLE VI: Calculations of Bhfs (in MHz) using new 
value and our Bhfa/Qnuc results reported in Table H] 



State 



B, 



hfs 



4p 'P3/2 
2 1 



5p 
3d 

4d 



"^3/2 



3/2 
3/2 



3d 2p'5/2 



4d 'Ds 



5/2 



-6.715(125) 
-2.260(56) 
-2.987(71) 
-0.774(16) 
-4.239(115) 
-1.100(24) 



fects in the present work. Combining our Bhfs/Q_ 



multi-configurational Dirac-Fock (MCDF) method [22 |. 
In contrast to our RCC approach, these methods account 
less correlation effects at the same level of excitations. 

As seen in Tables Hill and ITVl the trend of correlation ef- 
fects in A/j calculations in the first five low- lying st ates 
are same as discussed in our previous works [15l |35| ex- 
cept the differences in the magnitudes are due to the new 
basis functions. The amount of correlation effects in the 
higher excited states are comparatively smaller. Likewise 
in Ahfs, both the pair correlation and core-polarization 
effects which arise through Siy and RCC operators 
[15', 'ss'l, respectively, play major roles in obtaining the 
final results of S/i/s/Qrmc- The differences between our 
Bhfs/Qnuc result of the 3d state with others is 

due to the accurate treatment of these correlation ef- 



result with the measured B^fs value of this state 3lj . 
we get Qnuc — —0.0444(6)5. Following the same pro- 
cedure when we combine B^fs/Qnuc results of the 4p 
and 3d ^1)3/2 states with their corresponding ex- 
perimental Bhfs results, it gives less accurate values as 
Qnuc = -0.044(10)6 and Qnuc = -0.059(29)6, respec- 
tively. The associated larger errors are mainly due to 
larger uncertainties in the experimental results. We com- 
pare the above accurately estimated Qnuc value with the 
previously reported results in Table |Vl Our result is 
around 10% improvement over the considered standard 
value —0.049(5)6 [11] in this system. Recently Yu et al. 
[38t had evaluated this value as —0.0446, which is com- 
patible with our result, by combining their Bhfs/Qnuc 
result of the 4p state with its experimental Bhfs 

value which has around 20% uncertainty. Sundholm and 
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TABLE VIL Hyperfine interaction energies {Wf) for diflerent \F;I,J) states in ■*^Ca+ in MHz. 



F 


4s 51/2 


5s Si/2 


4p -Pl/2 


5p -Pl/2 


4p P3/2 


5p P3/2 


3d P3/2 


4d P3/2 


3d P5/2 


4d D5/2 


1 


















38.477 


32.834 


2 










201.400 


68.551 


317.553 


63.535 


32.444 


27.207 


3 


1814.371 


526.500 


327.201 


111.715 


115.087 


39.160 


177.841 


35.666 


22.940 


18.647 


4 


-1411.177 


-409.500 


-254.489 


-86.889 


-4.475 


-1.534 


-10.434 


-2.009 


9.421 


7.015 


5 










-161.121 


-54.824 


-248.977 


-49.932 


-8.840 


-7.880 


6 


















-32.752 


-26.271 



Olsen [2I] had combined precisely measured Bh fs of the 
3d 4s state of '^'^Ca with their calculated electric 
field gradient result to obtain Qnuc — —0.0408(8)6 in the 
same atom. In this work, they had employed the MCDF 
method on the restricted active space to calculate electric 
field gradient. We have already discussed the difference 
between the RCC and MCDF methods earlier in this sec- 
tion. There are also other works [13, [H, [H, [13, ES El 
finding Qnuc in '*'^Ca, but all of them have used either 
large uncertainty experimental results or more approxi- 
mated theoretical methods like second order many-body 
perturbation theory, results correcting for Sternheimer 
effects using the anti-shielding factor in the Hartree-Fock 
calculations etc. that cope with lesser electron correla- 
tion effects than our RCC method. 

Using our new Qnuc value, we determine B^fs from our 
Bhfs/Qnuc results presented in Table |T] and have given 
them in Table IVII The new Bh fs results are well within 
the error bars of experimental values with less uncertain- 
ties. Again, we evaluate energies of different hyperfine 
states using the formula given by Eq. (|2.2p correspond- 
ing to each fine structure level. In Table IVIIl we have 
reported these results which can be verified by analyzing 
isotope shift measurements in the future experiments in 
the assumed system. 



IV. CONCLUSION 

We have employed the relativistic coupled-cluster 
method in the Coulomb and Breit interaction approxi- 
mation to calculate the atomic wave functions in ^^Ca"*". 
Using these wave functions, we were able to determine 
Aiifs and Bhfs/Qnuc rcsults accurately. By combining 
our Bhfs/Qnuc result of the 3d ^£'5/2 state with its cor- 
responding precisely measured Bhfs value, we determine 
Qnuc of ^^Ca as —0.0444(6)5 which is 10% improvement 
over the considered standard value. In this work, we have 
also given contributions separately from the Breit inter- 
action and hyperfine interaction energies for a number of 
states. Our new Qnuc value and reported hyperfine inter- 
action energies may serve the researcher of both atomic 
and molecular physics. 

V. ACKNOWLEDGMENT 

This work is supported by NWO under VENI fellow- 
ship grant with project number 680-47-128. We are in- 
debted to B. P. Das for encouraging us to carry out this 
work. We also gratefully acknowledge discussions with 
C. Roos. The calculations were carried out using the 
Tera-flopp Super computer in C-DAC, Bangalore. 



[1] P. K. Ghosh, Ion Traps, Clarendon Press, Oxford (1995). 
[2] W. Demtroder, Laser Spectroscopy, Springer, Berlin 
(1998). 

[3] S. G. Karshenboim, S. I. Eidelmann, P. Fendel, V. G. 

Ivanov, N. N. Kolachevsky, V. A. Shelyuto, V. A.; T. W. 

Hansch, Nucl. Phys. B - Proc. suppl. 162, 260 (2006). 
[4] S. Boucard and P. Indelicato, Eur. Phys. D 8, 59 (2000). 
[5] C. Schwartz, Phys. Rev. 97, 380 (1955). 
[6] B. P. Das, K. V. P. Latha, B. K. Sahoo, C. Sur, R. K. 

Chausdhuri and D. Mukherjee, J. Theo. Comp. Chem. 4, 

1 (2005). 

[7] R. Kalish, Phys. Scr. 11, 190 (1974). 

[8] M. Keim, Recent Measurements of Nuclear Moments 
far from Stability, Exotic Nuclei and Atomic Masses 
(EN AM 98), International Conference at Bellaire, Michi- 
gan, USA, ed. by B. M. Sherrill, D. J. Morrissey and C. 
N. Davids, AIP Conference Proceedings 455 (Woodbury, 
New York), p. 50 (1998). 



[9] G. Neyens, Rep. Prog. Phys. 66, 633 (2003). 

[10] M. Arnold, J. Kowalski, T. Stehlin, F. Trager, G. Putlitz, 
Z. Phys. A 314, 303 (1983). 

[11] A. Wong, A. P. Howes, R. Dupree, M. E. Smith, Chem. 
Phys. Lett. 427, 201 (2006). 

[12] S. Forsen and B. Lindman, Methods of biochemical anal- 
ysis, ed. by D. Click, New York: John Wiley & Sons., 
27, 289 (1981). 

[13] C. Sur, K. V. Latha, B. K. Sahoo, R. K. Chaudhuri, B. 
P. Das, and D. Mukherjee, Phys. Rev. Lett. 96, 163003 
(2006). 

[14] B. K. Sahoo, Md. R. Islam, B. P. Das, R. K. Chaudhuri, 
and D. Mukherjee, Phys. Rev. A 74, 062504 (2006). 

[15] B. K. Sahoo, C. Sur, T. Beier, B. P. Das, R. K. Chaud- 
huri, and D. Mukherjee, Phys. Rev. A 75, 042504 (2007). 

[16] B. K. Sahoo, Phys. Rev. A 73, 062501 (2006). 

[17] B. K. Sahoo, Phys. Rev. A 74, 020501(R) (2006). 

[18] C. F. Roos, M. Chwalla, K. Kim, M. Riebe and R. Blatt, 



7 



Nature 443, 316 (2006). 
[19] H. Ito, Y. Li, M. Fujicda, M. Imac and M. Hosokawa, 
Optical frequency measurement precision of femtosecond 
laser optical comh system and the stability of its HF ref- 
erence frequency, Springer Series in Chemical Physics 79, 
846 (2005). 

[20] J. Benhelm, G. Kirchmair, C. F. Roos, and R. Blatt, 

Phys. Rev. A 77, 062306 (2008). 
[21] W. M. Itano, J. Res. Natl. Inst. Stand. Tochnol. 105, 829 

(2000). 

[22] W. M. Itano, Phys. Rev. A 73, 022510 (2006). 

[23] B. K. Sahoo, G. Gopakumar, R. K. Chaudhuri, B. P. 

Das, H. Merlitz, U. S. Mahapatra and D. Mukherjee, 

Phys. Rev. A 68, 040501 (R) (2003). 
[24] L. Wansbook, B. K. Sahoo, R. G. E. Timmermans, K. 

Jungman, B. P. Das and D. Mukherjee, Phys. Rev. A 

78, 050501 (R) (2008). 
[25] O. Lutz, A. Schwenk and A. Uhl, Z. Naturforsch. 28a, 

1534 (1973). 

[26] D. Sundholm and J. Olsen, J. Chem. Phys. 98, 7152 

(1993). 

[27] E. Bcrgmann, P. Bopp, Ch. Dorsch, J. Kowalski, F. 

Trager and G. zu Putlitz, Z. Phys. A 294, 319 (1980). 
[28] P. Grundevik, M. Gustavsson, I. Lindgren, G. Olsson, L. 

Robertsson, A. Rosen and S. Svanberg, Phys. Rev. Lett. 

42, 1528 (1979). 
[29] S. Salomonson, Z. Phys. A 316, 135 (1984). 
[30] R. E. Silverans, L. Vermeeren, R. Neugart, P. Licvcns and 

the ISOLDE Collaboration, Z. Phys. D 18, 351 (1991). 
[31] J. Benhelm, G. Kirchmair, U. Rapol, T. Korber, C. F. 

Roos, and R. Blatt, Phys. Rev. A 75, 032506 (2007); 

Phys. Rev. A 75, 049901(E) (2007). 
[32] Y. Ishikawa, R. Barctty and R. C. Binning, Jr. Chcm. 

Phys. Lett. 121, 130 (1985). 
[33] A. Mohanty and E. Clementi, Kinetically Balanced Ge- 
ometric Gaussian Basis Set Calculations For Relativis- 



tic Many-Electron Atoms, Modern Techniques in Cornpu- 

tatuinal Chemistry: MOTECC-89, ed. by E. Clementi, 

Chap. 4, pg. 169 (1989). 
[34] I. Lindgen and J. Morrison, Atomic Many-Body The- 
ory, ed. by G. Ecker, P. Lambropoulos and H. Walther, 

Springer- Verlag, Berlin (1985). 
[35] B. K. Sahoo, R. K. Chaudhuri, B. P. Das, S. Majumder, 

H. Merlitz, U. S. Mahapatra and D. Mukherjee, J. Phys. 

B 36, 1899 (2003). 
[36] B. K. Sahoo, G. Gopakumar, R. K. Chaudhuri, B. P. 

Das, U. S. Mahapatra and D. Mukherjee, Phys. Rev. A 

68, 040501(R) (2003). 
[37] B. K. Sahoo, Another paper will be reported shortly with 

an ab initio studies of electron affinity energies and po- 

larizabilities of many states using these basis functions. 
[38] K. -z. Yu, L. -j. Wu, B. -c. Gou and T. -y. Shi, Phys. Rev. 

A 70, 012506 (2004). 
[39] A. -M. Martensson and S. Salomonsson, Phys. Rev. A 

30, 712 (1984). 

[40] A. -M. Martensson, A. Ynnerman, H. Warston, L. Ver- 
meeren, R. E. Silverans, A. Klein, R. Neugart, C. Schulz 
and P. Lievens, Phys. Rev. A 45, 4675 (1992). 

[41] A. T. Coble and S. Maleki, Phys. Rev. A 42, 649 (1990). 

[42] F. M. Kelly, H. Kuhn and A. Pery, Proc. Phys. Soc. (Lon- 
don), Sect. A 67, 450 (1954). 

[43] F. Kurth, T. Gudjons, B. Hiblert, T. Reisinger, G. Werth 
and A. -M. Martensson, Z. Phys. D 34, 227 (1995). 

[44] W. Norterhauser, K. Blaum, K. Icker, P. Miiller, A. 
Schidmt, K. Wendt and B. Wiche, Z. Phys. D 2, 33 
(1988). 

[45] G. Breit, Phys. Rev. 34, 553 (1929). 

[46] P. Raghavan, At. Data Nucl. Data Tables 42, 189 (1989). 
[47] G. Olsson and S. Salomonson, Z. Phys. A 307, 99 (1982). 
[48] R. Aydin, W. Ertmer, U. Johann, Z. Phys. A 306, 1 
(1982). 



